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ABSTRACT 

The University of California Cerenkov detector on Pioneer 11 
previously observed Crand protons above 600 MeV in Saturn's inner 
magnetosphere, mixed with a poorly understood background of energetic 
electrons [Filllus and Mcllwaln, 1980]. Here we separate the electron from 
the proton counts and establish the first-order angular distributions for each 
species. To do this we use the theoretical relationships among the harmonic 
coefficients of the count rate as a function of spacecraft roll angle derived 
by Northrop [1985]. The majority of the counts were electrons with energy 
above several MeV; le, with drift periods shorter than the satellite orbital 
resonance. The electrons have isotropic pitch angle distributions, and the 
protons pancake over most of the region between Mimas and the rings, although 
there io a small region of dumbbell proton di’.stributlons in the vicinity of 
Janus and Epimetheus . 
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INTRODUCTION 

There wae an Internal Inconsistency In the data taken at Saturn by the 
University of California, San Diego Cerenkov detector on Pioneer 11. 

East-west anisotropies observed between the outer edge of the A ring at 2.28 
R and Mimas at 3.09 R apparently belonged to protons below the known 
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response range of the detector. In the original analysis, Fllllus and 
Mcllwaln deduced the gyroradlus by using Llouvllle's theorem to Interpret the 
radial displacement necessary to align the measured eastbound and westbound 
fluxes. The discrepancy was tentatively attributed to the presence of 
background electrons , but there was no method to quantify the contribution of 
these particles. Here we use the theory developed by Northrop In the 
accompanying paper to decompose the Cerenkov detector output Into two 
components, belonging to hlgh-rlgldlty protons and low-rlgldlty electrons. 

With this analytical tool we can give estimates of the proton and electron 
fluxes and, to a limited extent, their angular distributions. The results 
revise the high end of the energy spectrum of trapped protons created by the 
Grand mechanism, and reveal the existence of energetic electrons whose spatial 
distribution Is unlike that of the lower energy electrons. 

To show the Instrumental dilemma we first present an essential 
description of the detector and relevant calibrations . From the proton 
response we conclude that the east-west anisotropy was produced by protons of 
energy "600 MeV, Because the energy spectrum of the electron background is 
not well determined, we also Include a description of the electron response In 
order to evaluate two different possibilities , Then we review the 
observations briefly, and show how we apply Northrop"s method to separate the 
particle species. 
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PROTON RESPONSE OF THE CERENKOV DETECTOR 

The UCSD Cerenkov detector on Pioneers 10 and 11 has a radiator of 
60/40 methanol/water with index of refraction of 4/3. The liquid is held in a 
14 X 75 DiiD bottle of ultra-pure fused silica the front end of which is 
occupied by an expansion bellows for thermal compensation between the liquid 
and the bottle. The radiator is viewed at one end by a photomultiplier tube 
and associated electronics which count pulses above three discrimination 
levels (labeled Cl, C2, and C3) . No coincidence detectors are used because 
the high fluxes of Jupiter's radiation belts would have paralyzed the 
necessary circuitry. Instead, directionality is achieved by Internal 
reflection of light emitted toward the photomultiplier, and absorption of 
light traveling in the other direction on the bellows and blackened inner 
walls of the radiator housing. The discrimination levels are set so that the 
lowest channel (Cl) can be triggered by a particle whose pathlength in the 
radiator is not much more than a diameter, but the highest channel (C3) 
requires a pathlength comparable to the length of the radiator. Thus the 
channels have different angular responses, from Cl with significant side 
sensitivity, to C3 with a pencil lobe in the forward direction only. 

The angular response of the Cerenkov detector was demonstrated with an 
Identical unit in the proton beam of the NASA Space Radiation Effects 
Laboratory (SREL) synchrocyclotron. The detector was rotated and its 
effective cross-sectional area was normalized to that of a monitor detector. 
Figure 1 shows the effective area vs angle from the detector look axis for 
channel C3. 

When operating in a space radiation environment, the detector is 
irradiated simultaneously from all angles, and here the most useful 
calibration is the Integral of the effective area over solid angle, called the 
geometric factor, 

G(E) = / Agj^^(E) dn 

Figure 2 shows the geometric factor as a function of energy for several 
discrimination levels . Because the principal effect of varying the 
discrimination level is to vary the pathlength needed to yield enough light, 
the principal difference between channels is in the width of their angular 




responses, and so, In the magnitude of their geometric factors. The energy 
threshold is a second-order effect because all protons near the threshold 
easily have enough range to penetrate the entire detector and thus to produce 
the same pathlength. The channel to channel differences in energy threshold 
are primarily due to the fact that the lower channels are able to respond to 
Cerenkov light produced in the bottle walls and faceplate. For the index of 
fused silica, 3/2, the critical proton energy is 320 MeV. This accounts for 
the incipient response evident in channels Cl and C2 in Figure 2, but it is 
not a major factor because, with only a small amount of fused silica, it 
requir.'s exceptionally favorable particle trajectories. 

The peak energy of the SREL synchrocyclotron was not high enough to 
complete the response curve of the detector. To calibrate the asymptotic 
high-energy response, we used cosmic ray muon counting rates observed in our 
laboratory, corrected for the effect of the muon angular distribution over the 
angular response of the detector. These asymptotic geometric factors are 
plotted on the right-hand axis of Figure 2. 

The channels labeled Cl, C2, and C3 in figure 2 were matched to the 
prelaunch discrimination levels of the Pioneer 10 and 11 Cerenkov detectors. 
Both of these detectors underwent subsequent gain changes because of the high 
level radiation during the Jovian encounter. The effect on the Pioneer 
Cerenkov detector was to Increase the photomultiplier tube gain, which is 
equivalent to lowering the discrimination levels. The channel in Figure 2 
marked C3" is Interpolated from C2 and C3 to the level of Pioneer 11 channel 
C3 during the Saturn encounter. 

The small lump at "100 MeV is caused by scintillation of protons that 

stop in the optical materials. As mentioned by Fllllus and Mcllwaln J1980], 
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this represents an omnidirectional geometric factor of < .002 cm for 

protons between 65 and 140 MeV. By comparison, the UCSO solid state detector 
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on Pioneer 11 has a geometric factor of .012 cm over a larger energy window 
between 80 and several hundred MeV. Then by comparing the counting rates of 
these two sensors, we can conclude that no more than about 12%, of the C3 
response is caused by proton scintillation. 


ELECTRON RESPONSE OF THE CERENKOV DE-ECTOR 


The Cerenkov detector is also sensitive to electrons, (Its primary 
mission was to measure the Intensities of high energy electrons in the Jovian 
environment.) Figures 3 and 4 show the measured electron response. For 
electrons, unlike protons, penetration range is an important factor in 
understanding the detector response. The first Important threshold is the 
energy (“0.7 MeV) required for an electron to penetrate the detector housing. 
Since electrons are rela*;lvistic at this energy, the next criterion is the 
length of their track in the radiator, which determines their Cerenkov light 
output. As with protons, the particle trajectory must point toward the 
back end of the radiator in order to direct the Cerenkov beam to the cathode. 
The electrons" propensity to undergo high-angle scattering collisions has two 
effects. One is that the range tends to be understated by tables which list 
"end-point values," or the range of a scatterfree particle, and the other is 
that the angular responses tend to be smeared. 

Figure 3 demonstrates the net result for channel C3. Full response 
occurs at high energy where the particles" range exceeds the length of the 
radiator and the relatlvistli mass gain reduces their tendency to scatter. 

The biteout along the axis has two causes. One is that the asymptotic 
Cerenkov emission angle equals the critical reflection angle, so that light 
starts to be lost. The other is that there is more passive shielding at the 
front end of the housing than on the sides. As the electron energy decreases, 
the dominant response occurs on the side, where the radiator has a larger 
cross-sectional area, but smaller depth. Ultimately, as the electron energy 
decreases, the pulse height falls below the discrimination level. Then the 
electrons can only trigger the detector if several arrive simultaneously so 
that their pulses add in height. In this mode the response is peaked to the 
side where the cross-sectional area is largest. Channels Cl and C2 operated 
in this mode during the most Intense segment of the Saturn flyby. 

Figure 4 shows the geometric factor for single-particle events , 
plotted as a function of energy. This represents the detector's response to a 
delta function energy spectrum. It is convenient to represent this profile as 
a step function, so that we can quote the flux above some threshold. We 
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previously used Van Allen's *'bow~tie method" to obtain such a representation 
applicable to the Jovian radiation environment iFlllius and Mt’Ilwain, 1974], 
However, at Saturn, the electrons causing the pileup on channels Cl and C2 
have an extremely soft spectrum, outside the range of anything encountered 
previously. As the bowtie method sho'wed, the same threshold value is not 
applicable for all spectra. Therefore, .we have listed In Table I threshold 
energies applicable to a range of power law energy spectra, of the form, 
dN/dE " K e"" 

This table should help In interpreting the detector's response to very soft 
electron spectra. 


Table 1 

Threshold energies for channel C3 

usable for different power law energy spectra 

2 

with a geometric factor of 0,5 cm -sr 


n 

"th 

1 

14 

2.5 

14 

4 

12 

6.3 

10 

10 

8.5 


OBSERVATIONS 

The channel C3 counting rate reached a maximum in the vicinity of 2,7 
R (IR “ 60,000 km) and fell off to zero and near-zero values at the outer 
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edge ^f the A ring (2.28 R ) and the orbit of Mimas (3.09 R ). The east-west 
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anisotropies that we are Investigating are associated with the gradients on 
either side of the peak. Using a least-squares fitting procedure, Flllius and 
Mcllwaln [1980] represented the angular distribution of the counting rate by a 
truncated Fourier expansion of the foim C = Cq + CjSin(x) + C2Cos(2x) + 


C^elnCSx) + X being the angle between the look direction of the 

detector and the magnetic field i which happened to lie in the scan plane. 

Terms were chosen for symmetry and economy; omitted terms would violate mirror 
symmetry, and coefficients above 3x were so heavily convolved by the width 
of the sampling interval that they could not be evaluated with our limited 
amount of dots. In this representation the coefficient Cq Is the 
spln<-averaged count rate; the omnidirectional count rate is Cq - C 2/3 - 
Cj,yi5. The coefficient C 2 reflects whether the peak of the distribution is 
perpendicular to the magnetic field (C 2 < 0) or parallel to it (C 2 >0). The 
former case is called a pancake distribution and the latter, a dumbbell. The 
east-west anisotropy appears in the odd harmonics and C^. 

The coefficients derived from the data are shown in Figures 5 and 6 

for the Inbound and outbound radial cuts of Pioneer 11, The smooth curves 

were merely drawn by eye through the points. Note that outbound and 

have quite similar shapes, and Inbound they are similar, too, although the 

peaks are not quite coincident. In Figure 6 there was a data gap between 2,38 

Rg and 2.49 R^. Ihls has been filled in by interpolation so as to be able to 

carry out integrals over the entire range. Results inside 2.49 R outbound are 

8 

therefore not to be trusted and will not be given in later figures. 


SEPARATION OF PROTON FROM ELECTRON COUNTS 

The observed coefficients are assumed to be sums of proton and 
electron components. Thus C^ = C^^ + C^^, and similarly for the higher 
harmonic coefficients. We assume that there is no pile-up in this channel and 
no counts lost to dead time. In the preceding paper Northrop [1985] derived 
relationships among the harmonic coefficients of tvSe proton and electron 
count rate as a function of roll angle. Because the electrons have very small 
gyroradiuB , their odd harmonics vanish, and the odd harmonics observed are 
attributed entirely to protons. The even proton coefficients are given by 
equations (14) and (15) of that paper: 
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where is some point at which the count rate vanishes (say the outer edge of 
the A ring), and PpC^j) Is the proton gyroradius at r^ Based on the 
detector calibrations (See Figures 1 and 2), we used a proton energy of 600 
MeV to calculate the gyroradius . The Integrals were performed by Simpson's 
rule upon the smooth curves. This Is a bit subjective, but, unlike 
differentiation, integration tends to even out random errors in the input. 

Once Cqp and C2p are known, and are given by Cq^ = - 0^^ and 02^ = 


Coefficient is absent from Figures 5 and 6 because there were not 
enough data to determine it directly as a function of r. We have assumed that 
has the same shape as the other even coefficients and determined the ratio 
b « C^/Cq from the integral conditions (16) and (17) of the preceding paper. 
The conditions use the fact that the fluxes vanish at both boundaries of the 
region of interest. This method gives two values for b, one from each 
equation, and unless they agree, there is a dilemma. We find that inbound b ~ 
0.0915 and 0.0881 from (16) and (17) respectively. This minor discrepancy is 
handled by using 0.0915 in (1) and 0.0881 in (2) — that is, by using b in the 




equation from which It was determined. One may prove that the value of (1) or 
(2) is Independent of whether one integrates from T 2 Inwaru, or from rj^ 
outward to the r of Interest, provided that the b used in the equation was 
determined from It. Outbound the two values of b determined from (1) and 
(2)dlsagree: (1) gives b •• 0.00467 and (2) gives b ■ -0.0130. In this case we 
have again used the b In the equation from which it came. The disagreement 
may lie a consequence of the data gap. 

Figures 7 and 8 give the omnidirectional count rates for protons and 

electrons, and Figures 9 and 10 are the pitch angle distributions to 

lowest order in gyroradius. (The fg from equation 3 of the preceding paper.) 

We Immediately see that the proton count rate was only Z0% of the electron 

rate. The proton angular distribution Is found to be pancake (except In the 

vicinity of Janus and Eplmetheus where they arc dumbbell) and the electrons 

are Isotropic. Both angular distributions agree with what one expects from 

cosmic rays as the source, striking ring material in the case of the protons, 

or the moons Mimas, Janus, Eplmetheus, and rings for electrons [Blake et al., 

1983). That cosmic rays produce protons by the Grand process In the Inner 

region of Saturn’'s magnetosphere was first suggested by Filllus et al. [1980] 

and also by Cooper and Simpson [1980], Van Allen et al. [1980], and McDonald 

et al, [1980]. Other papers have subsequently studied the pitch angle 

distribution. If that distribution Is modeled by sln’"( j) , a value of 

n=6 at 3.5 R has been reported by Vogt et al, [1982]. Krlmigis and Armstrong 

[1982] find n=5, Schardt and McDonald [1983] find that for 48-160 MeV protons 

n increases from "1 at 2.7 R to nearly 5 at 2.85 R . 

s s 

DISCUSSION: INBOUND DATA 

The peak of the proton omnidirectional flux occurs at larger radial 
distance than the electron peak. (There is no theoretical reason for the 
peaks to coincide.) That there should be a proton peak Is to be expected; 
Grand is a distributed source of protons, and the A ring, Janus-Eplmetheus , 
and Mimas are sinks for the protons, which diffuse both ways from the Interior 
of the region under study. Cooper [1983] has solved the problem of the Grand 
source plus diffusion to absorbing edges and fit the University of Chicago 


data from Pioneer 11. From Figure 7 the peak for 600 MeV protons Is at 2.72 
R . and there Is a minimum at 2.52 R_ caused by proton absorption by Janus and 
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Eplmetheus, leaving a secondary peak at 2.4 R that Is about 1/3 the height of 

s 

the main peak* The position of the main peak and the factor of 1/3 both 
coincide with the profiles determined by Van Mien et al {1980] and Fllllus 
and Mcllwain [1980] for energies > 80 MeV. There Is no discernible absorption 
by the G-rlng at these >600 MeV energies, in contrast to lower energies where 
all experiments see an effect (See, for example. Van Allen [1983]). 

The proton angular distributions (Figure 9) are pancake outside the 
Janus-Eplmetheus absorption, and are dumbbell In the minimum, as would be 
expected from preferential wipeout of equatorial particles. The upswing at 
low count rates of these angular distributions Is probably not real, nor are 
negative values • There is much scatter in these deduced distributions , which 
Is not surprising considering the long train of analysis. The proton 
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distribution goes approximately as sin (6) at the peak (2.72 R ), but we 
cannot determine any trend at larger radii, although theory predicts more 
isotropy [Cooper, 1983]. 

In contrast to protons, there Is no significant volume source of 

electrons of the proper energy within the region. Two different sources of 

electrons are possible. High energy electrons are produced by cosmic rays at 

the rings and at Mimas, [Blake et al., 1983] and diffuse toward the center of 

the region from these sources at the two edges. The energy of these electrons 

is tens of MeV, and so they are well up on the plateau of the detector 

response (See Figure 4). The strength of this source Is unknown, and It Is 

not clear that there should be a peak at 2.62 R , given the location of the 

s 

sources • 


Also, electrons from other sources, outside Mimas, may be diffusing 

inward. The spectrum of these electrons is soft, because the preferred energy 

for electrons to escape absorption by Mimas Is just 1 MeV, and by 7 MeV the 

relative drift time has fallen from infinity to <5 hours. As compiled by 

Chenette and Stone [1983], the electron spectrum at Mimas falls by 4 orders of 

magnitude between 5 MeV and 11 MeV. Magnetic-moment-conserving diffusion of 

relativistic electrons from Mimas to the electron peak at 2.62 R_ Increases 
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the energy of electrons with 90° pitch angle by a factor of 1.3, With the 
steeply falling energy spectrum i this energization would increase the number 
of electrons above the energy threshold by a factor of 21. The count rates at 
Mimas were small , at most a few per second » so it seems that a peak of 313 per 
second is coo large Co be accounted for completely in this way. 

In sum, we do not know what fraction of the electrons are produced by 
cosmic rays at Mimas and Che rings, and what fraction by other processes 
elsewhere. 


OUTBOUND DATA 

The electron peak is comparable in height with the inbound electron 
peak but is at a larger radius than J Abound , As Inbound, the electron angular 
distributions are quite Isotropic. Tlie radial profile of the electrons is 
different from that observed at lower energies, in that it is peaked, and 
vanishes at the orbits of Mimas and at the F ring. At lower energies (“1 MeV) 
the electrons do not exhibit clear channels at the orbits of these moons, 
which is consistent with their resonant energy to avoid sweeping. 

The peak in proton flux in Figure 8 is at the same radius (2,72 R ) as 

O 

Inbound and is at just about the same height. As at lower energies (e.g., 

above 80 MeV (Van Allen et. al, 1980; Fllllus and Mcllwaln, I960]) there 

appears to be azimuthal symmetry about Saturn. Although we have no data 

inside 2.49 R , there is no evidence in Figure 8 of an incipient upturn of the 
6 

curve Inward of 2,5 R^, Van Allen does see wings of a secondary peak outbound 
that he estimates would have been lower than Inbound by a factor of 2 to 3. 

He suggests that such an inbound-outbound asymmetry may be connected with the 
longitude of moons such as Janus and Eplmetheus relative to that of Pioneer 
Inbound and outbound. 

The proton angular distributions are pancake except near Janus and 
Eplmetheus where they are dumbbell — just as was the case inbound. The 
pancake distributions are less well organized than Inbound. We do not know 
whether this is real or merely represents the noise in the- analysis, which is 


not completely eatlsfylng h$ic.ausc of the determlnisttlon of C4 and the data 
gap# NcvcrthclcBB , the resultB of the anelyslB ogrec remarkably well with the 
Inbound pasB( and with obaervotlons ot lower energies (>80 McV). 


CONCLUSIONS 

!• We have employed the method developed by Northrop In the 
accompanying paper [19B5]> for deducing the omnldlrcctlonol Intensity and 
first-order angular distributions of hlgh-rlgldlty particles In the presence 
of a background of low-rlgidlty particles. The method Is applied by 
performing appropriate Integrals over the observed angular distributions! and 
can be thought of as an Integration of the radial gradient deduced from 
the east-west anisotropy of the hlgh-rlgldity particles. 

2. We have evaluated the number of 600 MeV protons necessary to 
produce the cast-west anisotropy observed by the UCSD Cerenkov detector on 
Pioneer 11 between 2.3 and 3.1 R^. This Is a significant downward revision of 
the estimate made by Filllus and Mcllwaln [1980 1', but nevertheless still 
consistent with other data at lower energies. Crand is the most likely source 
for these particles, as concluded by Filllus and Mcllwain [1980]. 

3. ’There exists a high energy component of the electron spectrum 
which has a spatial profile with low values at Mimas and Janus -Eplmcthcus and 
a maximum in between. This profile indicates that sweeping takes place at 
these moons, which is to be expected for electrons of several MeV or more. 

The source and spectrum of these electrons remains uncertain. 

4. We have determined the pitch angle distributions of these 
particles to first order in gyroradlus. The electrons are almost isotropic, 
but the protons exhibit a combination of pancakes and dumbbells, with the 
latter near Janus and Eplmetheus . 
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FIGURE CAPTIONS 


Figure 1 

Effective Area of Cerenkov detector channel C3 as a function of angle 
to the detector axis for protons of energy 66, 98, 137, A46, 497, 521, and 560 
MeV. The proton energy Is given at the right hand edge of each curve. 

Because of experimental difficulties there Is significant scatter and 
statistical error, and this Is illustrated for 560 MeV (circles) and 98 MeV 
(crosses). These data were taken at the NASA Space Radiation Effects 
Laboratory ( SREL) , 

Figure 2 

Response vs energy for several pulse height channels of the UCSD 
Cerenkov detector. (For Identification of the channels, see the text.) The 
proton response was obtained at SREL up to its peak energy of 560 MeV, and 
ground level cosmic ray muons were used to calibrate the asymptotic response 
of the detector to particles with ■> 1. 

Figure 3 

Effective area vs angle for electrons of the Indicated energies.. 

These data were obtained on the electron llnac operated by Intelcom Radtek 
Corporation In San Diego. 

Figure 4 

Response vs energy for channels Cl, C2, and C3 of t’»e UCSD Cerenkov 
detector. To represent these profiles as Ideal step functions, one must allow 
the value of the threshold energy to depend upon the energy spectrum of the 
incident radiation. (See Table I and text.) 

Figure 5 

Harmonic coefficients of the count rate as a function of the 
spacecraft roll angle, plotted vs. radial distance to Saturn during the 
Inbound pass . 

(a) Even harmonics, Inbound. 

(b) Odd harmonics, Inbound. 




Figure 6 

Harmonic coefficicnta of the count rate as a function of the 
spacecraft roll angle • plotted ve. radial distance to Saturn during the 
outbound pass. 

(a) Even harmonics, outbound* 

(b) Odd harmonics, Inbound. 

Figure 7 

Derived omnidirectional count rates for protons and electrons on the 
Inbound pass . 

Figure 8 

Derived omnidirectional count rates for protons and electrons on the 
outbound pass . 

Figure 9 

Derived angular distributions for protons and electrons on the inbound 
pass , 

Figure 10 

Derived angular distributions for protons and electrons on the 
outbound pass • 
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